from pylab import *

nbin = 20

w0 = -1.0
wa = 0.0

zmax = 1.5
amin = 1./(1+zmax)

#z = linspace(0,zmax,20)

a = linspace(1.,amin,20)
z = 1/a-1

w = zeros(z.shape)

for i in range(len(z)):
	w[i] = w0 + wa*z[i]/(1.+z[i])

#print 'DDE_z = ', 
#for i in range(len(z)):
#	print '%5.3f, '%(z[i]),


#print ''

#print 'DDE_w = ', 
#for i in range(len(w)):
#	print '%5.3f, '%(w[i]),

for i in range(len(w)):
    print 'DDE_w'+str(i) + '= %5.3f %5.3f %5.3f'%(w[i],-2.5,0.5)
    
print around(z,5)
